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O ' Abstract 

00 : 

O ■ We establish new convergence results, in strong topologies, for so- 

lutions of the parabolic-parabolic Keller-Segel system in the plane, to 
the corresponding solutions of the parabolic-elliptic model, as a phys- 
ical parameter goes to zero. Our main tools are suitable space-time 

5^ \ estimates, implying the global existence of slowly decaying (in gen- 

eral, nonintegrable) solutions for these models, under a natural small- 
ness assumption. 

1 Introduction 

We consider two related nonlinear parabolic systems which are frequently 
used as models for a description of chemotactic phenomena, including the 
aggregation of microorganisms caused by a chemoattractant, i.e. a chem- 
ical whose concentration gradient governs the oriented movement of those 
microorganisms. The parabolic character of the systems comes from the 
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diffusion described by tlie Laplacians. A version of tlie system (PE) be- 
low is also used in astrophysics as a model of the evolution of a cloud of 
self-gravitating particles in the mean field approximation. 

The first one is the classical parabolic-elliptic Keller-Segel system 



X e M^ t > 0, (PE) 



Here, u = u{x,t), ip = (p{x,t) are either functions or suitable (tempered) 
distributions. When u > 0, ip > 0, they may be interpreted as concentra- 
tions (densities) of microorganisms and chemicals, respectively. 
The second one is the parabolic-parabolic system 

Ut = An — V ■ (uVip), 

Tipt = Aip + u, a;GM^t>0, (PP) 

u{0)=uo, <^(0) = 0, 

where r > is a fixed parameter. Each of the models can be considered 
as a single nonlinear parabolic equation for u with a nonlocal (either in x 
or in {x,t)) nonlinearity since the term Vif can be expressed as a linear 
integral operator acting on u. In the latter model, the variations of the 
concentration ip are governed by the linear nonhomogeneous heat equation, 
and therefore are slower than in the former system, where the response of (p 
to the variations of u is instantaneous, and described by the integral operator 
(— A)^^ whose kernel has a singularity. Thus, one may expect that the 
evolution described by (PE) might be faster than that for (PP), especially 
for large values of r when the diffusion of ip is rather slow compared to that 
of u. Moreover, the nonlinear effects for (PE) should manifest themselves 
faster than for (PP). 

The theory of the system (PE) is relatively well developed, in particular 
when this is studied in a bounded domain in W^, d = 1, 2, 3, with the homo- 
geneous Neumann conditions for u and (p at the boundary of the domain. 
One of the most intriguing properties of (PE) considered for positive and 
integrable solutions u in d = 2 case is the existence of a threshold value Svr 
of mass M = J u{x, t) dx, see the pioneering work [17J and [HE]. Namely, if 
Mo > is such that / Uo{x) dx > Svr, then any regular, positive solution u of 
(PE) cannot be global in time. We refer the reader for a fine description of 
the asymptotic behaviour of integrable solutions of (PE) in the subcritical 
case M < Stt to [H] and to ^U\ for the limit case M = Svr. See also [8], [9] 



in the radially symmetric case. The higher dimensional versions of (PE) 
have been also extensively studied, cf., e.g., [21 El [5], and [1] for blow up 
phenomena. 

The doubly parabolic system (PP) has been a bit less studied. For 
instance, it is known that if for the initial data Uq one has M < Svr, then 
positive solutions are global in time, see [12] in the case of a bounded planar 
domain, and also [3] Theorem 5]. However, it is not known whether M < 
8tt is, in general, a necessary condition for the existence of global in time 
solutions. That is, it is not known whether the blow up occurs for solutions, 
except for a specific example in [16] of a particular blowing up solution for 
a system close to (PP). Even, it is an open question what is the exact range 
of M guaranteeing the existence of integrable self-similar solutions. For the 
system (PE) it is proved that M G [0,87r), and the self-similar solutions 
(unique for a given M G [0, Stt)) describe the generic asymptotic behaviour 
of global in time, positive and integrable solutions of (PE). Concerning (PP), 
it is known that M < M{t) with M{t) linear in r, is a necessary condition 
for the existence of self-similar solutions, cf. [1]. For (PP) with small M 
such special solutions are also important in the study of space-time decay of 
general solutions, see [2l|. The analysis if any M > Svr may correspond to a 
self-similar solution is under way, see [7|. For a different point of view about 
self-similar solutions for higher dimensional models of (PP), see also |18j . 
Usual proofs of a blow up for (PE) involve calculations of moments of a 
solution and then symmetrization, cf. [TH [T]. These methods seem do not 
work for (PP), hence another approach is needed to show a blow up for that 
system. For a numerical insight on blow-up issues we refer, e.g., to |14j . 

A nice result in [25j shows that the solutions of the systems (PP) and 
(PE) enjoy a kind of stability property as r \ 0: solutions of (PP) converge 
in a suitable sense to those of (PE). It had been an old question raised by 
J. J. L. Velazquez and D. Wrzosek, recently solved in [23]. However, this 
result obtained for suitably small solutions in quite a big functional space of 
pseudomeasures, gives no indication on the behaviour of possible ( "large" ) 
blowing up solutions. 

The solvability of the systems (PE), (PP) has been studied in various 
classes of functions and distributions, like Lebesgue, Morrey, Besov, etc., 
with an immediate motivation to include the a priori strongest possible 
critical singularities of either solutions or initial data which appear to be 
point measures in the two-dimensional case and the multiples of \x\~'^ func- 
tion in the higher dimensional case. In particular, "vast" functional spaces 



suitable for analysis of the two-dimensional systems include measure and 
pseudomeasure spaces, cf. plISllGj 125]. 

We show in this paper a result on the existence of (in general, nonin- 
tegrable) solutions in a class X of functions with natural space-time decay 
properties, see Theorem 12.11 and 13.11 Here, the space E of admissible initial 
conditions also contains Dirac measures. The corresponding solutions may 
be positive and "large" in the sense of their nonintegrability. Nevertheless, 
they are defined globally in time. Unlike the paper P3], we work in {x,t) 
space, when [25] has dealt with the Fourier variables ^, cf. the formulation 
( l36l) below. Such results are obtained by an extension and refinement of 
techniques used in [S] for (PE) in higher dimensions, but neither for (PP), 
nor in the two-dimensional case of (PE) which often requires a specific 
treatement. Moreover, the function spaces that we employ here allow us to 
deal with data that can be more singular than those considered in [5] . The 
spaces X and E defined in the next section are, in a sense, critical for that 
analysis, and have been already considered, in slightly different forms, e.g. 
in the studies of the Navier-Stokes system in [131 ES] . 

Our main results are contained in Section |U where we address the prob- 
lem of the convergence as r \ of solutions m^ of the system (PP) to the 
corresponding solutions u of (PE), in the space X, arising from small data 
in E. Mathematically, our stability result is not included in, and does not 
imply, that of [25]. However, it seems to us that the use of the natural (x, t) 
variables provides a more immediate physical interpretation. Furthermore, 
our method looks more flexible, and can be used to prove the stability of 
the system with respect to stronger topologies. For example, we establish 
also the convergence in the Lf^(L^)-norm for data belonging to EnL^, and 
in the L^-norm for data in ECiL'^. Motivated by [19], we will also address 
this issue in the more general setting of shift invariant spaces of local mea- 
sures. The main difficulty for obtaining the convergence u"^ ^ u in strong 
norms is that Vip'^ enjoys some kind of instability as r — > 0, in particular 
in weighted spaces. 

Moreover, we give a nonexistence (blow up) result for solutions of (PP) in 
W^, d > 1, with the positive Fourier transform of Uq in the spirit of [23], see 
Theorem 15. 1[ These are complex valued solutions with no straightforward 
physical/biological interpretation. However, such a result tells us that there 
is no hope to prove the global existence of solutions to (PP) and similar 
models for arbitrarily large data relying only on size estimates. 



2 The parabolic-elliptic system 

In order to study the systems (PE) and (PP) we introduce the Banach space 
X of functions u = u{x, t) and the Banach space E of tempered distributions 
Mo G 5'(M^) by defining the norms 

ll'"IU = esssupt>o,^gM2(t+ \x\'^)\u{x,t)\, (1) 

and 

\Wo\\e = lle'^^MolU- (2) 

Here, e*^ denotes the heat semigroup defined by the Gaussian kernel gt, 
gt{x) = 4^e~4tl^l . For example, the Dirac mass Uq = 5 in M? is an element 
of E. Notice that, by the definition, E is continuously embedded into the 
weak Hardy space 7i^, which is the space consisting of all tempered dis- 
tributions / such that sup(>o \e^^f\ belongs to the Lorentz space L^'°°(]R^). 
See [21] • 

Let us define the bilinear form Bq by 

Bo{u, v){t) = [ e(*-^)^V ■ {uV{-A)-^v){s) ds. (3) 

Jo 

Here, (— A)^^ is the convolution operator on functions defined on M^ with 
the kernel K{x) = — ^log|x|. With this notation, the equivalent integral 
(mild) formulation to (PE), called also the Duhamel formula, reads 

u{t) = e^^uo- Bo{u,u). (4) 

We begin by establishing the following simple result 

Theorem 2.1 There exist two absolute constants e, f3 > with the following 
property. Let Uq & E be such that \\uo\\e < e- Then there exists a unique 
(mild) solution u E X of (PE) such that \\u\\x < e/3. 

The proof of Theorem 12.11 will follow from a series of lemmata. 
First , we have the following estimate of the leading term in Vy? 
Lemma 2.1 Let u E X and yj be such that Aip + u = 0. Then 

Vifi{x,t) = ^ f uiy,t)dy + nix,t) (5) 

Fl J\y\<\x\/2 

with Co = — 2~ (I'lT'd the remainder TZ satisfying 

\n{x,t)\ < c\\u\\x (t^ + \x\ 



Proof. Indeed, let us represent the partial derivatives of ip, for j = 1, 2, as 

djip = —-J *u = Ii + I2 + h, 
where 

J\y\<\x\/2 \x-y\^ 

The terms I2 and Is are obtained by taking the integration domains {|a; — 
y\ < |2;|/2} and {|x — y\ > \x\/2, \y\ > \x\/2}, respectively, in the convolu- 
tion integrals defining dj(f. It is straightforward to prove that I2 and J3 can 

be bounded by C||n||;f (t^ + |x| j . On the other hand, we can rewrite /i 



^ f uiy,t)dy + R^ix,t). 



as 

Co: 

kP J\y\<\x\/2 

An application of the Taylor formula shows that the above bound holds also 
for Ri. □ 

We immediately deduce from (^ the following useful estimate 
Lemma 2.2 Let u E X and A(f + u = 0. Then 

||V<^(t)||L- <C\\u\\xt--^. 

The last lemma that we need is the following 
Lemma 2.3 Let u,v E X . Then, for some constant Co independent ofu,v 

\\Bo{u,v)\\x < Co\\u\\x\\v\\x- 

Proof. We can assume, without any restriction, that ||'u||a' = II'^IIa' = 1- 
Lemma 12.21 implies 

\uVi-A)-^v\{x,t)<C\x\-h-\ (6a) 

and also 

\uV{-A)-^v\{x,t) <C\x\-h-i. (6b) 

We denote the gradient of the heat semigroup kernel gt by 

G(..<).V.(^e-iHA. p) 



Then we may represent Bq as 

Bo{u,v){x,t)= f fGix-y,t~s)iu\/i-A)-'v)iy,s)dyds 
= Ji + J2, 
where Ji = /o4|<|^|/2 • • • and J2 = IoI\y\>\^\/2 ■■■ ■ Using the estimate 

\Gix-y,t-s)\<C\x-y\-lit-s)--^ 
and inequahty (lUal) . we get the bound 

\Ji{x,t)\<C\x\-'^. (8) 

Another possible estimate is 

\Ji{x,t)\<C\x\-h-^, (9) 

which is obtained using the bound 

\G{x-y,t-s)\ <C\x-y\-'^{t-s)-^. 
On the other hand, from the property 

||G(-,t-s)|Ui=c(t-s)-5 (10) 

and inequahty ( l6bl) . we obtain 

\J2{x,t)\<C\x\-\ (11) 

As before, we have also the bound 

\J2ix,t)\<C\x\-h-^. (12) 



This second estimate is deduced from (!6al) . 

Then, using ((H]), flTTl) . we obtain the space decay estimate 

\Boiu,v)\{x,t)<C\x\~^ (13) 

and from ([9]), ( TT2l) — a provisory (not optimal) time decay estimate 

\\Bo{u,v)it)\\j,^<Ct--K 

But we may represent Bq as 

Boiu,v){t) = e^^Boiu,v){t/2) + f G{t - s) * {u\/i~A)-^v){s) ds. 



Jt/2 
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Thus, applying the weak Young-type inequahty for convolutions in Lorentz 
spaces L3'°°*L'*'^ C L°°, see [19], and the equality obtained from the scaling 
laws in Lorentz spaces, 

1 



AlTt 



J_| |2 

-e «' 



_3 
Ct 4. 



(14) 



L4.1 



we finally get 

\\Bo{u,v){t)\\L^ < Ct-^ + C f {t-s)-^\uV{-A)-^v{s)\\L^ds 

Jt/2 



(15) 



Combining inequalities (fT3l) and (fT5l) we get Bo{u, v) G A", together with its 
continuity with respect to u and v. □ 

Proof. Note that (using the duality S — S') we have c^^uq -^ uq in S' as 
t — >■ 0. The conclusion of Theorem 12. II follows in a standard way (cf., e.g., 
P^ [21 E]) from the contraction fixed point theorem. □ 



3 The parabolic-parabolic system 

Let r > be a fixed parameter. We consider the system (PP) whose 
equivalent integral formulation reads 



u{t) 



e'^uo 



V • e(*-^)^ 



u 



^ ^ Jo 



ds. 



(16) 



We introduce for all r > the bilinear form B^ (recall that G is defined by 
the expression ([7])) 

B^{u,v){x,t) = f f Gix-y,t~s){uWr{v)){y,s)dyds, (17) 



where Wr{v) is the linear operator acting on v 



Wriv)ix,t) 



T 



G{*-^]*v{a) 



(x, a) da for r > 0, 



:i8a) 



with a natural convention 

Wo{v){x,t) = {V{-A)-'v) {x,t). (18b) 

In this way, the system (PP) is also rewritten in a compact form (cf. (jlj)) 



as 



u = e*^Uo~Br{u,u). (19) 

We are going to solve flT^ in the space X exactly as was in the parabolic- 
elliptic case. An additional estimate, however, is needed: 
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Lemma 3.1 Let u E X and r > 0. Then there exists a constant C* > 0, 
independent of u and r, such that 



\WM{t)\\L^ <CH-nu\\x. 



(20) 



Proof. As usual we can and do assume \\u\\x = 1- Then, for all 1 < p < cxd 
we have 



— — — i-i-i 
\u{x, a)\ < \x\ pcr p . 



This implies, for 1 < p < oo. 



I / Ml ^ -1 + - 



Now, we represent 



where h = J^ ... and I2 = j\ 



\h 



L°° 



<c 



. Evidently, we obtain the bound 

t/2 



t 

t/2 



W-' 



T \ T 



I ■u(cr) 11^2,00 da 



L2.1 



<ct-^. 

For the integral I2, let us begin with a rough bound 



21'-; L° 



<C 



t/2 



i^n-a 



T \ T 



Infer) ILoo dcr 



(21) 



< Cr^^ri 



This bound gives the required estimate, excepted when r belongs to a neigh- 
bourhood of the origin. Thus, in the sequel, it is enough to consider the 
case < r < i. Now, we further decompose 



I2 = hi + I: 



2,2, 



where /2,i = f^,^^ ■ ■ ■ and 12,2 = Jt-rt • • • • Next, we are going to improve 
(J2T]) writing 



1^2,1 



<c 



t-Tt 



t/2 



Ig"-" 



T \ T 



t-Tt 



7 1 



H(j)\\j,^da 



< Cre / (t - a) 6(7 3 da 

Jt/2 

< ct-^. 



(22) 



The last estimate is 



\h,2{t)\\L^<C I 

Jt-Tt 



i_^n-a 



T \ T 



\\u{a)\\Loc da 



< 



Cr-5 f {t-a)-^a-^da 

Jt-rt 



(23) 



<cri 



The conclusion of Lemma O follows from (EH), (122]), (ESD- □ 

Lemma Em allows us to see that, ii u, v & X with \\u\\x = \\v\\x = 1, 
then 

\uWr{v)\{x,t) <C\x\-h--i, (24a) 

and 

\uWr{v)\{x,t) <C\x\-H-^, (24b) 

for some constant C > independent of r > 0. These are the analogous 
estimates as those we obtained in the parabolic-elliptic case (see inequal- 
ities ( !Ua|) and (l6bl) ). Then, using exactly the same arguments as in the 
previous section, we arrive at the following existence result 

Theorem 3.1 There exist two absolute constants e*,l3* > with the fol- 
lowing property. Let uq ^ E be such that \\uo\\e < e*. Then there exists 
a unique (mild) solution u & X of (PP) such that \\u\\x < t*(3* . 

Remark. The case of nonzero initial data v^(0) can be studied in a quite 
similar way. 

Remark. A closer look at the proofs of estimates for Vy? in (PE) and 
(PP) reveals that the behaviour of Vy? is a bit different in these two cases. 
Namely, if ^ -u and u{x,t) ~ (t+ ja^P)""^ (in the sense that Ci(t+ Ix]"^)^^ < 
u{x,t) < C2(t + la^n^^ for some Ci, C2 > 0), then it follows from ([5]) 
that Vip{x,t) ~ la;!^"*^ log(t + jxp) for (PE), while Vip is more regular: 

|V(y5(a;, t)| < c (t2 + |a;| j in (PP) case. In other words, letting y be the 
space of functions / = f{x, t) such that /^ G X, we have (p'^ E y ior r > 0, 
but {(f'^) does not converge in 3^ as r — > 0. However, such an instability does 
not prevent from the convergence of the densities u'^ ^ u for vanishing r. 
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4 Study of the r \ limit 

We now study the convergence as r \ of solutions u'^ of the system (PP) 
to the corresponding solution u of (PE). A result in this direction has been 
obtained recently by A. Raczyhski in [22], who established the convergence 
u^ ^ u in the norm 3^^, for a G (1, 2), defined as 



esssup (l+t^|e|)"|M(e,t)|. (25) 

We will obtain in subsection 14.21 a similar result using the A'-norm. 



t>0,,eGK2 



4.1 Regularity properties of solutions of (PE) 

In this subsection we prepare some preliminary material. The first Propo- 
sition consists of a regularity result with respect to the space variable for 
solutions of (PE). The second Proposition describes their regularity prop- 
erties with respect to the time variable. 

Proposition 4.1 For all r G (1, 2) there exists a constant e,., with < e^ < 
e (the absolute constant of Theorem \2.1\) such that, if \\uq\\e < ^r, then the 
solution of (PE) constructed in Theorem \2.1\ satisfies 

||VM(t)||^,,^<Crt+^ (26) 

for some constant C = C{uo,r) independent oft. 

Proof. We use a standard argument involving the subspace A^ C A" defined 
by 

Xr = lue X, 3C: \\Vu{t)\\Lr,^ < crt+H, 

and equipped with its natural norm. Recalling that the kernel f^ of the op- 
erator V(— A)^^ belongs to L^'°°, first we deduce from the Young inequality 

\\V'{-ArMt)L..^ < Ct-'^^'\\v\U, i = i - i (27) 

Next, from the Holder inequality (noticing that ||M(t)|| 2,2,00 < Ct^2 HmH^^^,), 

\\uV\-A)-'v{t)\\Lr.oo < Ct-'+^\\uy\\v\\;,^. (28) 

The generalization of the classical inequalities to Lorentz spaces can be 
found, e.g., in 
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We claim that the bihnear operator Bq introduced in (^ is boundedly 
defined: Bq: Xr x Xr ^ Xr. Indeed, for \\u\\xr = \\'v\\xr = 1; we combine 
the estimates 

||VG'(t-s)||i<C(t-s)-i 

and (TTUl) with the inequahty (a consequence of m G A") 

\\u{s)\\Lr,o. <Cs-^+^, (29) 

the estimate || V(— A)^^f (s)||oo < Cs^2 obtained from Lemma l^^ and 
Then we arrive at 

\\VBo{u,v)(t)\\Lr:oo 

pt/2 

<C / {t - s)-'\\u{s)\\Lr.o.\\V{-/\)-^v{s)\Uds 
Jo 






3 + 1 



Moreover, for uq G E, we have \e2^uo{x)\ < C(t + Ixp)^"*^. Hence, 
||e2^Mo||L7,oo < Ct~ ^^ for 1 < 7 < oo. We now choose /3,7 G (l,C)o) such 
that 1 + ^ = ^ + -. Then, the semigroup property of the heat kernel gt, 
and the fact that Vgt/2 G L^'^, imply 



Now the usual the application of the contraction mapping theorem, in 
a closed ball of small radius in the space Xr, allows us to conclude. □ 

The following proposition is the first crucial tool for our stability re- 
sult. It provides the Holder regularity, with respect to the time variable, of 
solutions of (PE) in Lorentz spaces. 

Proposition 4.2 Let 1 < r < 2 and Uq G E, such that \\uq\\e < ^r- Then 



the solution u of (PE) constructed in Proposition 4J_ satisfies for all < 
t' < t 



1 , .,. _3 , 1 



\\u{t)-u{t')\\Lr,oo <C{t-ty{tT^^~, (30) 

for some C = C{uo, r) independent oft and t' . 

Proof. It is enough to show that both c^^Uq — e* ^Uq and BQ{u,u){t) — 
BQ{u,u){t') satisfy the required bound in the L^'°°-norm. 
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From the identity 
e^^uoix) - e^'^uoix) = / \e^'^uo{x ~ y) - e^''^uo{x))gt-t'{y) dy 
= - j j Ve*'^Mo(x - Oy) ■ ygt-Ay) dydO, 
we get 

Now, we can write 

Bo{u,u){t) - Boiu,u){t') = Ai + ^2, 

with 

Ai= f {G{t -s)- Git' - s)) * (mV(-A)-^m) (s) ds 
Jo 



A2= f G{t - s) * (mV(-A)-^m) (s) ds. 
Jt' 



and 



Recall that from m G ^Y we deduce fl29|) . Combining this with the estimate 
of Lemma 12.21 we get 

\\uV{-/\)-'u{s)\\^^^^<Gs-i+-r. (31) 

This immediately yields 

IIA2IU..00 < f (t-s)-^s-t+^ds<C(t-t')^(i')~^^^- 
Jt' 

The estimate of Ai is slightly more involved. We start with the identity 

G{t -s)- G{t' -s) = ('e^*-*')^ - Id) G{t' - s). 

The action of the convolution operator with the function on the right hand 
side is studied via the following variant of a result established in [20] . 

Lemma 4.1 Denote the Calderon operator by A = (— A)2 . Then, for some 
constant C depending only on r G (1, 00) 



e*^-Id)/||^.,^<Ct^||A/lU,.,^. 
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Proof. Writing / = A ^A/, we see that 

(e*^-Id)/ = $,*(A/), 
where, 

and 

$(o = (e-i«|'-i)ier^ 

It only remains to show that $ G L^(R^), which is immediate. Indeed, it 
is well known, and easy to check, that ^{^) = |^|e~'^' defines a function 
\E' G L^(]R^) (for example, with the method described in [15], one obtains 
|*(x)| < C(l + |x|)~^ and |V^(x)| < C(l + |x|)~^). We conclude applying 
the Bochner inequality to the identity 



/oo 
\E'(?7x)d?7. 



Using this Lemma we deduce 

||Ai||i.,oo < C{t - t')^ / ||AV^t/_, * {uV{-Ay^u) {s)\\Lr.oo ds 

Jo 

= C(t-t')^(A,l+A,2), 

where Ai^i and Ai 2 are obtained splitting the integral at s = t'/2. 

But, as the function \E' introduced in the proof of Lemma 14.11 satisfies 
V^E' G Z/^(M^), we see by a simple rescaling that 

\\AVgt'-s\\i<C{t'-sy\ 
Combining this estimate with inequality ( !3T|) . we get 
Ai,i<C{t')-^^K l<r<2. 
For treating Ai^2, we combine the estimate 

\\Agt'^s\\i<C{t'-s)--^ 
with the inequality (for 1 < r < 2) 

ll(Viz) (V(-A)-Im) (s)|U..oc + \\u {\/\-A)-'u) (s)|U..oo < Cs-^-^'r, 

obtained by applying ( l26l) . Lemma [2^21 and ( l28l) with u = v. We get as 
before 

Ai,2<Cit')-'i+-r, l<r<2, 

and this concludes the proof of Proposition 14. 2[ □ 
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4.2 The vanishing r Umit 



After Proposition I4.2[ the second crucial step for the study of the hmit as 
r \ consists in the asymptotic analysis of the linear operators Wr, t > 0, 
introduced in fll8ap - fll8b|) . This is the purpose of the following lemma. 

Lemma 4.2 Let e = e{T) be an arbitrary function, strictly increasing and 
continuous on [0, 1], such that e{0) = 0. Let also 1 < r < 2 and u be 
a function satisfying, for < t' < t. 



(32) 

\\u{t)-u{t')\\Lr,o.<C{t-t')-2{t')--^''-r, 

with a constant C independent oft,t'. Then, for allt>0,T& [0, 1], and 
for another constant C , independent of t and t, 

sir) supti \\(WAu) - Woiu)) (t)|L < Cr^-i (33) 

In particular, if uq E E is small enough (for example, H^^olU ^ ^3/2^ then 



the corresponding solution u of (PE) constructed in Proposition\4.1\ satisfies 



hm supt^ \\{WAu) - Woiu)) (t)|L = 0. (34) 

Proof. Without any restriction we can assume that < £:(r) < | for 
positive T. Define e(r) such that e = e~^2 . Borrowing from [25] the idea of 
splitting the time integral using intervals depending on r, we write 



where 



next 



and 



Wr{u) - Wo{u) = Ji + J2 + J3 



ds. 



Ji{t) = j -G[-y-)*u{s) 

J2{t) = [ -gC-^^) *u{t)ds - Wo{u){t), 

Jt(l-e(r)) T \ T J 



Mt) = [ -gC—^) * [u{s) - u{t)] ds. 

Jt(i-e{r))r V r J 

From the first relation of (1321) and the Young inequality in Lorentz space 

(using that, by ([7]), G{-,t) G L'' '^(M^), where r' is the conjugate exponent), 

we get 

/•t(l-e(r)) 

||Ji(t)||oo<Cr^-^ / {t - s)-"^-ls-^+-r ds 

Jo 
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Notice that this estimate of Ji is exactly what we need for 

As for J2, we see from a simple computation via the Fourier transform 
that 

J2{t) = -9te(r)/r*Wo{u){t). 

If - = - — ^, then we deduce from the usual weak-convolution estimates that 
Wo{u) = V(— A)~^u is bounded in the L°'°°-norm, by Ct~^^~. Applying 
once more the Young inequality (using now gte(T)/T ^ L"' '^), we get, as 
before, 

\\j2mo0 < c i^j^y' %-i 

Applying the second of inequalities ( !32l) . we obtain immediately 

\\Mt)\\oo<CT'--k-K 



which is even better than what we need. This proves the inequality (133|) . 
Choosing, for example, r = | and £(r) = |t^^^^ proves the last claim (jSll) 
of Lemma 14.21 □ 

We are now in the position of establishing our first main result 

Theorem 4.1 There exists an absolute constant e' > (a priori smaller 
than the constants e, e* > in Theorems \2.1\ and \3.1\) . such that if Uq G E, 
II Wo lis < e'; then denoting by u E X the solution of (PE) and u^ E X the 
solution of (PP) constructed in the previous theorems, we have as r \ 

u^ ^f u in X. 

Proof. The proof follows easily from Lemma 14.21 Indeed, from the integral 
equations (jlj) and (IT9l) . the bilinearity of Br and i?o, and the smallness of 
the solutions u^ and u, we have (similarly as in [25j|, where two terms in the 
bihnear expansion can be absorbed by the left hand side) 

Wu" -u\\x < C\\Br{u,u) - Bq{u,u)\\x. 

But, by the definition of Wr and Wq (see (ITSaD-lfTSbl)). 

Br{u,u){t) - BQ{u,u){t) = f G{t-s)*{u{Wr{u)-Woiu))){s)ds. 

Jo 



Argueing as in the proof of Lemma 12.31 we obtain 



\Br{u,u) - Bo{u,u)y < C\\u\\xhnpt-^\\Wr{u) - Wo{u){t)l 
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If e' > is small enough, then Lemma 14.21 can be applied to the solution u 
of (PE), implying that the right hand side of the above inequality has a 
vanishing limit for small r. This finally gives 

ll^"^ — mIIa:" ^ as r \ 0. 



Remark 4.1 Notice that, the smaller the norm ||tio|U; the faster the con- 
vergence u'^ -^ u as r — > 0. This is due to the fact that for very small data 
it is possible to apply Lemma 112] with r close to 1 (despite the constants in 
our estimates blow up as r \ 1). More precisely, our arguments show that 
for any < 5 < |, one can find a constant C > and e{6) > such that, 
for ||mo||e < e((5), one has Hu"^ — u\\x < Cr^^^ for all < r < 1. 

4.3 The parabolic-elliptic limit in stronger topologies 

If Mo > is small in the E'-norm, and belongs to a smaller space, for 
example, Uq G E H L^, then the solutions u'^ and u of (PP) and (PE) will 
remain in L^, uniformly in time, during their evolution. Hence, it is natural 
to ask whether the convergence u^ ^ u holds also in the natural norm 
of L°°((0, oo); L^). Our next theorem provides a positive answer. As the 
proof of this fact does not really depend on a particular topology under 
consideration, it seems appropriate to consider a more abstract setting. 

We denote by C any shift invariant Banach space of local measures, see 
p^ Ch. 4] for their definition and main properties. These are Banach 
spaces of distributions, continuously embedded in P'(R^). Moreover, they 
are known to satisfy the following properties (for some constant C > 
depending only on £), 

1. For all / G £, g E L^(M^), the convolution product f*g is well defined 

in£and||/*(7||^<C||/|Ulklli. 

2. For all / G £, /i G L°°{E?), the pointwise product fh is well defined 

in/:and||//i|U<C||/|U| 



3. Each bounded sequence {fk} C C has a subsequence convergent in £, 
in the distributional sense. 

Obvious examples of spaces satisfying these properties (and which are in- 
deed shift invariant space of local measures) are the L^-spaces, 1 < p < oo, 
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the Lorentz spaces L^'"^, l<p<oo, l<g<oo and the space of bounded 
Borel measures M(M^) = Co(M^)*. In the latter case, such duahty rela- 
tions ensures Property 3. Other interesting examples include the Morrey- 
Campanato spaces M^''', (1 < p < g < oo) and suitable multiplier spaces, 
see UHl Ch. 17]. 

On the other hand, the space of pseudomeasures, i.e., the space of tem- 
pered distributions / such that / G L°°) does not fulfill the second re- 
quirement. Therefore, the stability result in the pseudomeasure topology 
will not be encompassed by our next Theorem, but requires a specific (and 
more involved) treatment, see [2^ . 

Notice that, because of the conservation of the total mass for positive 
solutions of (PE) and (PP), the L^-norm remains constant during the evolu- 
tion. This observation will allow us to handle the case of data uq & E (1 L^ , 
despite Property 3 breaks down for L^. 

Theorem 4.2 Let uq E E r\ C, where C is either a shift invariant Banach 
space of local measures, or C = L^ . In the latter case we require either 
Uq > or, for signed Uq, that ||mo||i is sufficiently small. Then there exists 
a positive constant I , depending only on C, such that, if 

\\uq\\e < e, 

then the systems (PE) and (PP) possess unique solutions u and u^ , respec- 
tively, such that for an absolute constant (3 > 0, \\u\\x < /3e and \\u'^\\x ^ /3e. 
In addition, 

sup ||M(t)||£ < oo and sup ||M'^(t)||£ < oo. 

t>0 i>0 



Moreover, the conclusion \\u'^(t)—u(t)\\x -^ of Theorem\4.1\is strengthened 



sup \\u^ it) — u{t)\\c ^0 as t\0. 
t>o 

Proof. Obviously we have, for some constant Cq > independent of t 

We'^uoWc < Co. 

Moreover, for each r > (we include in this way the analysis of (PE)), we 
have the estimate 

\\B^iu,v){t)\\c < C(snp \Ht)\\c) Wvy, (35) 



where C > depends only on C. This follows from flT7|) written as 

Br{u,v){t)= f G{t-s)*{uWr{v)){s)ds 
Jo 



for each r > 0, with the convention (118bl) . Therefore, using the usual 
estimate ( ITOl) . we have 



\\B^{u,v){t)\\c<Cisnp\\u{t)\\c](snpt-^\\WAv){t)l 

^t>0 ^ ^t>0 

The last factor is bounded by \\v\\x owing to Lemma [2.21 in the case r = 0, 
and to Lemma [2.31 for r > 0. This yields ( l35l) . 

Now, we can consider, for r > 0, the sequence of approximating solutions 

ul = e*^Mo - Br{ul_^, ul_^), k = l,2,... . 

When e < min{e, e*}, we know by the proofs of Theorems 12.11 and 13.11 that 
the sequence u^ converges in X to the solution u^ of (PE) or (PP). Here, of 
course, u = vP for the solutions of (PE). 

On the other hand, applying recursively ( l35l) . we get ulit) G C for all k 
and 

^^V\\ul{y)\\c < Co + pCeisnp \\ul_^{y)\\c) , 

with P = max{/5, /5*} (the constants obtained in Theorems 12.11 and 13. Ih . 
Iterating this inequality we arrive at 

snp \\ul{t)\\c <C' <oo, 

provided /5Ce < 1. If £ is a shift invariant Banach space of local measures, 
from Property 3 we get for all r > 

sup ||M'^(t)||£ < C" < oo, 
t>0 

where C" > is independent on r. Of course, the last claim remains valid 
in the case C = L^ and Uq > (notice that the smallness of ll^olU prevents 
blow up results that could occur, otherwise, when the second moment of mq 
are finite and j Uq > Svr. See the introduction and the references therein 
quoted). If we remove the assumption uq > 0, we can obtain the same 
conclusion provided ||mo||i is sufficiently small. Indeed, we see from inequal- 
ity (IH31) that the fixed point argument applies in the space Lf{L^) fl X. 
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We now discuss the stability, including also the case C = L}. From the 
bilinearity of Br-, the mixed estimate (l35l) and the smallness of the solutions 
u and u^ (allowing two terms of the bilinear expansions to be absorbed by 
the left hand side), we obtain 

\\u^{t)-u{t)\\c<Csn^\\Br{u,u)-Bo{u,u){t)\\c. 

t>0 

Argueing as at the end of the proof of Theorem 14.11 we arrive at 



\\Br{u,u)-Bo{u,u){t)\\c < C(snp\\u{t)\\c] isnpr^\\{Wr{u)-Wo{u)){t)l 
and the conclusion follows from Lemma HI 



Remark. As an application of this general result, let us observe that taking 
C = L°°, we obtain for Uq & E Cl L°°, with Uq small in the i?-norm, 

as r \ 0, uniformly in (x, t) G M^ x [0, oo). 

5 Blow up of complex valued solutions of the 
parabolic-parabolic system 

Consider the system (PP) in the space M."^ with any d > 1. Passing to the 
Fourier transforms, we may write the Duhamel formula (iT6l) in the form 

miUle'^t-sMl^^^^is^aM^^^ - r], s)u{r], a) dr] da ds. 

(36) 

Our goal is to construct a class of complex valued initial data, such 
that the corresponding solutions blow up in finite time, in any classical 
norm. For a G M, we denote by B^ the homogeneous Besov space, which 
can also be identified with the Holder-Zygmund space C". As it is well 
known (see fi9[ l23]). most of the classical functional spaces (including all 
homogeneous Triebel-Lizorkin, and thus Lebesgue or Sobolev spaces) are 
continuously embedded in -B^°° for some real a. 

Theorem 5.1 There exists wq G 5(M'^) such that woiO — ^> ll'^olUi = 1- 
// A is sufficiently large, then any (tempered) distributional solution of (1361) 
with Uq = Awq (and thus any solution of (PP)) blows up in a finite time. 
More precisely, for some time t* < oo and each a ^ M., \\u{t*)\\^a,oo = cxd 
holds. 
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Our approach is closely related to that in ^ Theorem 3.1] which followed 
the argument in [23] for the "cheap" Navier-Stokes equations. We produce 
some estimates from below of the Fourier transform of any solution with 
u{0) = uq that can be obtained via the iteration procedure for (l36ll with 



Jo ^l Jo 



for k = 1,2,... . This recurrence relation, in general, does not preserve 
the positivity of the Fourier transform of Uq. This leads us to restrict 
our attention to data of the form uq = Awq with wq G 5(]R'^), such that 
supp WqC {^eR'^: 2-1 < 6 < l-^l < !}• Define, for A; = 0, 1, . . . the set 

Then we see that for w^ = Wq , we have Wk = {2TT)''^Wk-i * Wk^i, and 
therefore, suppw/; C Ek- This implies that if, in addition, wo{^) > 0, then 
the positivity of the Fourier transform will be preserved by the sequence Uk, 
and so by u. Next Lemma tells us more: 

Lemma 5.1 For all k = 0,1,2, .. . , we have 

u{i,t) > [3ke-^"'lt,<t<At)wk{i), (37) 

where {(3k} and {tk} are two sequences defined below in (l38l) and (1391) . 

Proof. For A; = 0, the conclusion immediately follows from 

provided we choose (3q = A and to = 0. Let A; > 1. Assume that the 
inequality of the lemma holds for A; — 1. Then, for all tk < t < t* , 

t s 

Wk-i{^ - V)wk-i{v) . . . 
X — -T dri da as 

(271)^^ ' 

t 



tk-i 



Wk-ii^-v)wk-M . . 
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Thus, we can bound u{C,, t) from below as follows 

t 

92^-3 



u[ 



{U)> / is-t,_^)- e-(*-^)2^'=e-^(*'-*-^)2^'=-VLie-^^-^^.(0d. 



t 

This suggests us to set, for some 5 > 0, t* — t^^i = 5r2"^'''+^. With this 
choice, putting to = 0, we have 

t* = 5T, tk = 5T{l-2-^^). (38) 

Then, tk - tk-i = 35t2-'^''. We get, for tk < t < t*, 

u{^,t) > {36t - 1 + e-''n'^-''-'T-'e~'pl_,e-'''wk{0- 

This inequality is interesting only if the right hand side is positive. Therefore 
we will assume that 36t > 1. We choose {/3fc} in such a way that 

/5o = A, Pk = {36t - 1 + e-^^^)2-2^-V-ie-^/3ti, k = l,2,... . 

This choice leads to inequality (1371) . 

In order to compute Pk, we introduce Ms^r, such that 

2^^- = {36t - 1 + e-^^^)e-^ 2-\~\ 

Notice that we have jSk = 2^^^'^~'^^Pl_-^^. We claim that, for some a,b,cE M, 

/5fc = A2'2"+''^+'=2'=_ 

Indeed, from an easy calculation we find b = 2, a = A — Ms,r and finally 
c = Ms^r — 4, which is needed to ensure Pq = A. 
Hence, we have 

Pk = (A2^''--^y' 2^-*^^.-+2fc^ A; = 0, 1, . . . . (39) 
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We conclude that when 

we have by flHUl) /3fc ^ oo and, in particular, ||u(tfc)||/,i = ||M(tfc)||oo -^ oo for 
k ^ oo. The above size condition on A can be rewritten in an equivalent 
form as 

A further analysis of the lower bounds obtained for the Fourier transform of 
a candidate solution permits us to conclude, as was in [S], that ||M(t*)||^a,oo = 
cxD for each a G M, so that all Besov (and also L^ or Triebel-Lizorkin) norms 
of u blow up not later than t*. Notice that, for a blow-up at t* = 1, we need 
A > Ce^/^r, cf. 
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